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We study various methods of constructing fat links based upon the HYP (by Hasenfratz & Knechtli) and Fat7 
(by W. Lee) algorithms. We present the minimum plaquette distribution for these fat links. This enables us to 
determine which algorithm is most effective at reducing the spread of plaquette values - a strong indicator of 
improved statistics for spectrum and static potential measurements, among other quantities. 



1. INTRODUCTION 

Despite its prominent role as a useful tool 
in the non-perturbative study of QCD, lattice 
field theory has long struggled with the difficul- 
ties of accurately representing a continuum field 
theory in a discrete spacetime. The occurrence 
of discretization errors has lead to the develop- 
ment of a great many improvement schemes, some 
more useful than others. One extremely popular 
method of improvement has been the use of "fat 
links" , in which each ordinary or "thin" link of 
the lattice, U^{x), is replaced by a linear com- 
binaton of the thin link and the links adjacent 
to it. Two particular fattening algorithms will 
be of interest to us in this report, HYP or hy- 
percubic blocking £Q, and Fat7 blocking (that is, 
SU(3)-projected Fat7 blocking 0). These two 
algorithms have been shown to be perturbatively 
equivalent at one-loop level |2I3| . We will present 
the results of a preliminary numerical investiga- 
tion to determine whether there are any signifi- 
ant non-perturbative differences between the two 
algorithms. These results will point the way to 
more detailed investigations. 



2. HYPERCUBIC BLOCKING 

The HYP algorithm was proposed as a mod- 
ification of standard APE smearing 0j. APE 
smearing involves constructing a fattened link 
V fl (x) by adding the sum of staples, weighted by 



some factor a, to the thin link, 

V^x) = (l-a)U fi (x) + ^J2 W ^ W 

followed by projection back to SU(3) (where 
W v (x) refers to the sum of three-link staples 
in the positive and negative v directions, with 
central link parallel to U^x)). Clearly n re- 
peated applications of such an algorithm will 
access the gauge fields at a distance n lattice 
spacings from the original thin link (although 
smoothing algorithms such as cooling and smear- 
ing can be viewed as diffusive "random- walk" pro- 
cesses which primarily affect physics only within 
a range a^/n, as investigated in 5 ). To lo- 
calize APE smearing within the smallest possi- 
ble volume, HYP blocking was defined as being 
equivalent to three iterations of APE smearing 
with the caveat that the staples at each stage 
could not be constructed from links which are 
fattened in the same plane as the staple itself. 
The original thin link is therefore fattened by 
pairs of positively-and- negatively-oriented staples 
in three directions. Each of these staples is con- 
structed from links fattened by staples in only 
two (positive and negative) directions, and each 
of these staples is constructed from links fattened 
in a single plane (for further details please re- 
fer to £Q). HYP blocking therefore uses only the 
staples defined by the boundaries of hypercubes 
attached to the original thin link (Fig^). Projec- 
tion back to SU(3) is performed after each level 
of fattening. 
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Figure 1. HYP fattening. The bold line is the 
thin link, and the dashed lines are the hypercubes 
attached to it. The double and zig-zig lines are 
the staples at the final and middle levels of fat- 
tening respectively. The outer level (fourth direc- 
tion, which must necessarily be constructed first) 
is omitted. Only positively directed staples are 
shown here. 



We shall henceforth refer to the links fattened 
by the outermost set of staples as the initial-stage 
fat links, the links fattened by staples formed 
from the initial-stage links as the middle-stage fat 
links, and the completely fattened links (i.e. those 
fattened by staples formed from middle-stage fat 
links) as the final-stage links. 

Since its introduction, HYP blocking has 
shown itself to be extremely effective at reducing 
taste-symmetry breaking effects, and has gained 
widespread attention and use. 



3. FAT7 AND FAT7 BLOCKING 

Fat7 blocking (i is a form of fattening which 
incorporates not only standard three-link staples 
but also five- link and seven-link staples (Fig. 0. 
Clearly, with staples up to length seven it is pos- 
sible to traverse all three directions orthogonal to 
the original thin link, as is the case with HYP 
blocking. However Fat7 blocking as it was orig- 
inally envisaged did not incorporate SU(3) pro- 
jection. Fat7 is a modification of this algorithm 
which incorporates SU(3) projection of the final 
sum of link and staples ■ 

If we construct Fat 7 fat-links in an analogous 
manner to the iterative construction of HYP links 
we have greater freedom in our application of 
SU(3)-projection. An iterative approach allows 



Figure 2. The thin link plus 3-link, 5-link, and 
7-link staples used to construct the Fat7 fat link. 
In practice the thin link and each of the staples 
would be weighted by an appropriate prefactor. 



us to include SU(3) projection at each level of 
the construction of the link, rather than just at 
the final level. The only difference between HYP 
and Fat7 is then that in HYP blocking each link 
in each staple is fattened at the subsequent level, 
whereas in Fat7 it is only the links parallel to the 
original thin link that are fattened (see Fig. |2J. 

We can see that the second covariant derivative 
operator may be used to recursively define a 
fattened link 

L v {a) ■ M x ) = (1 - 2 ") ' U„(x) + aW u {x) (2) 

(W v (x) defined as in Eq. Q). With the param- 
eter a taking the value one-quarter this oper- 
ator can be interpreted as suppressing flavour- 
changing gluon interactions by vanishing in the 
limit as gluon momentum approaches the lattice 
cut-off. The Fat 7 link can easily be constructed 
by repeated application of this operator as 

V » = I E ■ (L P (a) ■ (L x (a) ■ UJ) (3) 

pcrm(^,p,A) 

where we note explicitly that the sum over permu- 
tations of the directions is taken. This repeated 
use of the same operator accounts for the use 
of a single weighting parameter a at each stage, 
whereas HYP blocking uses three parameters a\, 

a 2 , a 3 . 

As noted above, to construct the Fat7 links 
from the Fat7 links, SU(3) projection must be 
performed after the final summation is taken, but 
it can also optionally be applied after each or 



Figure 3. The HYP fattening (left) compared 
with Fat7 (right). The original (thin) link is the 
bold black line, while the double, zig-zig, and 
solid thin lines represent the staples at each suc- 
cessive level of fattening. Only positive directions 
are depicted here, for the sake of clarity. 



either of the first two fattening stages, leaving 
us with four possible combinations of projection 
schemes to define our Fat7 link (i.e. projection 
at the final, initial+final, middlc+final, and ini- 
tial+middle+final stages). 

4. COMPARISON OF FATTENING AP- 
PROACHES 

Fat links have shown themselves to be ex- 
tremely useful at overcoming problems associ- 
ated with the lattice formulation. Specifically 
they have shown a reduction in the severity of 
taste-symmetry breaking associated with stag- 
gered fermions, and a reduction in the severity 
of the exceptional configuration problem associ- 
ated with the Wilson fermion formulation. Both 
of these problems can be directly related to short- 
scale fluctuations of the link values. In pQ, the 
values of the minimum plaquettes calculated on 
an ensemble of configurations were used as a 
probe of the most severe link fluctuations. The 
maximum plaquette values are bounded above, 
and hence we expect an increase in the minimum 
plaquette value of any configuration to indicate a 
smaller range of link fluctuations across that con- 
figuration, and hence a reduction in the scale of 
the largest short-range link fluctuations. 

The process of fattening tends to increase the 
minimum plaquette values of different configura- 
tions to differing degrees, spreading the minimum 
values over quite a wide range (as can be seen 
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Figure 4. Histogram of minimum plaquette values 
for our ensemble of unblocked configurations. 



by comparing Figs. 0] and Figs. and EJ). We 
expect that certain measurements will achieve a 
lower level of statistical uncertainty if performed 
on an ensemble of configurations having a narrow 
spread of minimum plaquette values. One of the 
aims of this preliminary study will be to identify 
whether the Fat7 blocking algorithm is capable 
of producing a narrower spread of minimum pla- 
quette values than HYP, and is hence worthy of 
a more detailed investigation. 

5. RESULTS 

We present the minimum plaquette values for 
an ensemble of 194 quenched configurations gen- 
erated on an 8 3 x 32 lattice at j3 = 5.7. We also 
present the minimum plaquette values for the en- 
sembles obtained after blocking each configura- 
tion with the HYP algorithm (Figs.[51and[f)|), and 
with the Fat7 algorithm (Figs-dtoHUJ). As noted 
above the natural choice for the Fat7 weighting 
parameter is a — 0.25. For HYP we are interested 
in two choices of parameters, those found non- 
perturbatively in pp, namely a\ = 0.75, ct2 = 0.6, 
«3 = 0.3 (HYP-I), and those found perturba- 
tively in 3 , namely a\ = 7/8, a.2 — 4/7, 0*3 = 
1/4 (HYP-II), although we note the caveat that 
choices of «i larger than 0.75 may tend to destabi- 
lize the smearing algorithm, making the blocked 
configuration rougher rather than smoother 8 . 
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Figure 5. Minimum plaquette values for HYP-I 
blocked configurations. 
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Figure 6. Minimum plaquette values for HYP-II 
blocked configurations. 
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Figure 7. Minimum plaquette values for Fat 7- 
blocked configurations, with final stage SU(3) 
projection. 



For the Fat7-blocked configurations we have anal- 
ysed the four different SU(3) projection options 
described in section |21 The histograms of min- 
imum plaquette values are presented in Figs. 01 
to [IT]] All figures are to the same scale. 

We can see clearly that the initial unblocked 
data produce a very narrow peak. After HYP 
blocking the distribution of minimum plaque- 
tte values has been shifted towards (but not 
into) the positive region, and taken on a Poisson 
distribution-like outline. The data obtained from 
Fat7 with SU(3) projection at the final stage, and 
at all three stages are especially interesting from 
an analytical basis 0. It is clear that all four 
Fat7 fattening schemes lead to minimum plaque- 
tte distributions that have a substantially nar- 
rower spread of values than the HYP-I and HYP- 
II data. Although the number of configurations 
used in this study is small, the results seem sig- 
nificant enough to warrant further investigations 
with Fat7. 



6. CONCLUSIONS 

Although this work is only preliminary, we find 
encouraging signs that Fat 7 may be capable of 
producing results which are as good as HYP, and 
possibly better for some calculations. Both HYP 
and Fat7 algorithms can be made substantially 
faster and more efficient by pre-calculating and 
storing the staples across the entire lattice, and 
using these to construct the links at the next level 
of fattening. This approach avoids calculating 
the same staple more than once, as one would 
be forced to do if one fattened in a naive site- 
by-site manner. Using this approach Fat7 can be 
made less memory-intensive than HYP, as it only 
requires the staples in a limited number of direc- 
tions to be pre-calculated at each stage of fatten- 
ing. An efficiently-designed algorithm for Fat7 
does not sacrifice speed, and enables us to use 
only one-quarter of the memory that is needed to 
store the pre-calculated staples for HYP. In future 
work we hope to determine whether this cheaper 
form of fattening merely equals, or in fact can 
produce results which are superior to HYP for 
certain calculations. 
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Figure 8. Minimum plaquette values for Fat 7- 
blocked configurations, with middle and final 
stage SU(3) projection. 
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Figure 9. Minimum plaquette values for Fat 7- 
blocked configurations, with initial and final stage 
SU(3) projection. 
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Figure 10. Minimum plaquette values for Fat 7- 
blocked configurations, with initial, middle, and 
final stage SU(3) projection. 



